###4.2. Cervejas

cervejas <-c("chope","lata","garrafa","chope","garrafa", "garrafa","lata",
"lata","nenhuma","lata","garrafa","garrafa", "garrafa","lata","lata",
"lata","garrafa","lata","chope","nenhuma", "garrafa","garrafa","garrafa",
"chope","garrafa","garrafa","chope","garrafa","lata","lata")


##4.2.1.Represente este resultado como um grfico de barras e um dotplot 
#(funo dotchart).

(tab.cervejas<-table(cervejas))
barplot(tab.cervejas)
dotchart(tab.cervejas)

##4.2.2.Qual tem maior razo dado/tinta?
#O dotplot tem a maior razo dado/tinta

###############

###4.3 Caixetais

##4.3.1.Construa um histograma do dap[,5] dos fustes dos caixetais

hist(caixeta[,5],xlab="DAP dos fustes",main="Histograma")
##se eu quissesse fazer por fuste
##unique(caixeta$fuste)
##hist(caixeta[caixeta$fuste==1,5])

##4.3.2.Construa histogramas da altura das rvores para os diferentes 
#caixetais ('local')

hist(caixeta[caixeta$local=="chauas",6],xlab="Altura",
main="Histograma Chauas")

hist(caixeta[caixeta$local=="jureia",6],xlab="Altura",
main="Histograma Jureia")

hist(caixeta[caixeta$local=="retiro",6],xlab="Altura",
main="Histograma Retiro")

#summary(caixeta[caixeta$local=="retiro",6])
#summary(caixeta[caixeta$local=="chauas",6])
#summary(caixeta[caixeta$local=="jureia",6])
min(caixeta[caixeta$local=="jureia",6])
max(caixeta[caixeta$local=="jureia",6])


##4.3.3.H diferenas entre as estruturas (distribuio de tamanhos) dos 
#caixetais?

#Sim, "chauas" tem uma maior distribuio de tamanho do que os outros 
#locais.  Enquanto "Retiro" apresenta a menor distribuio de tamanhos.

###############

###4.4. Eucaliptos
#utiliza o arquivo e.grandis
e.grandis<-read.table("egrandis.csv", header=T,sep=";")
str(e.grandis)
head(e.grandis)
summary(e.grandis)

##4.4.1.Utilize o grfico boxplot para analisar o DAP de rvores de E. grandis 
#em funo das variveis regio (regiao) e rotao (rotacao). 

boxplot(e.grandis$dap~e.grandis$regiao*e.grandis$rotacao)

##4.4.2. Avalie a normalidade da altura do conjunto total de rvores com um 
#grfico quantil-quantil contra a distribuio normal.

qqnorm(e.grandis$ht)
qqline(e.grandis$ht)

#A altura do conjunto total de rvores no apresenta distribuio normal,
#dados das pontas so diferentes do esperado pelo calculo terico#

###############

###4.5. Mais caixetais

###calcule as reas basais de cada fuste, e ento some-as
###deduza desta rea basal total o dap total
####altura media dos fustes

##4.5.1.Analise a relao dap-altura ('dap' e 'h') em funo do caixetal
# (local) com a funo plot, mas somente para as rvores de caixeta
# (Tabebuia cassinoides). 

caixeta<-read.table("caixeta.csv", header=T, sep=",", as.is=TRUE)
head(caixeta)
#xlab="Altura", ylab="DAP", main=tit.chauas
head(caixeta)
tab.caixeta<-caixeta[caixeta$especie=="Tabebuia cassinoides",]
head(tab.caixeta)
tab.caixeta$area.b.por.fuste<-pi*((tab.caixeta$cap/2)^2)
names(tab.caixeta)

area<-aggregate(tab.caixeta[,8],by=list
(local=tab.caixeta[,1],arvore=tab.caixeta[,3]),FUN=sum)
names(area)
area$dap<-sqrt((area$x/pi))
head(area)

altura<-aggregate(tab.caixeta[,6],by=list
(local=tab.caixeta[,1],arvore=tab.caixeta[,3]),FUN=mean)
names(altura)
altura$dap<-area$dap
names(altura)

#plot(area$dap~altura$x,xlab="Altura", ylab="DAP", main=tit.chauas)

tit.chauas<-"T.cassinoides em Chauas"
tit.jur<-"T.cassinoides em Jureia"
tit.ret<-"T.cassinoides em Retiro"

plot(area$dap~altura$x,subset=area$local=="chauas",
xlab="Altura", ylab="DAP", main=tit.chauas)

plot(area$dap~altura$x,subset=area$local=="jureia",
xlab="Altura", ylab="DAP", main=tit.jur)

plot(area$dap~altura$x,subset=area$local=="retiro",
xlab="Altura", ylab="DAP", main=tit.ret)

## 4.5.2.Para a mesma relao do item anterior, verifique linearidade com a 
#funo scatter.smooth

scatter.smooth(area$dap[area$local=="chauas"]~
altura$x[altura$local=="chauas"], col="red",
xlab="Altura", ylab="DAP",main=tit.chauas)

scatter.smooth(area$dap[area$local=="jureia"]~
altura$x[altura$local=="jureia"], col="blue",
xlab="Altura", ylab="DAP", main=tit.jur)

scatter.smooth(area$dap[area$local=="retiro"]~
altura$x[altura$local=="retiro"], col="green",
xlab="Altura", ylab="DAP", main=tit.ret)

## 4.5.3.Utilizando o pacote lattice, analise a relao dap-altura ('dap' e 'h') em funo do caixetal (local),
# mas somente para as rvores de caixeta (Tabebuia cassinoides).

library(lattice)
xyplot(dap~x|local, data=altura, ylab="DAP", xlab="Altura")
